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Abstract. "Streamers, sprites, leaders, lightning: from micro- to macroscales" was 
the theme of a workshop in October 2007 in Leiden, The Netherlands; it brought 
researchers from plasma physics, electrical engineering and industry, geophysics 
and space physics, computational science and nonlinear dynamics together around 
the common topic of generation, structure and products of streamer-like electric 
breakdown. The present cluster issue collects relevant articles within this area; most 
of them were presented during the workshop. We here briefly discuss the research 
questions and very shortly review the papers in the cluster issue, and we also refer to 
a few recent papers in other journals. 
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1. Streamers: questions and methods 

1.1. The pivotal role of streamers in discharges in matter 

Rapid electrical breakdown is a generic and ubiquitous process in nature and technology; 
its earliest stage is governed by the formation of filamentary streamers. They occur 
in lightning in the earth's atmosphere and on other planets, they are developed for 
and used in various technical and industrial applications, and they are studied as a 
fundamental phenomenon in air, pure nitrogen, argon, and other gasses. They occur 
equally in water, oil, semiconductors and other materials. Streamers play a pivotal role 
in the complex sequence of events that start with electron avalanche and culminate in 
leader or arc discharges, and their internal processes are the most nonthermal within this 
sequence. While the various stages of electrical breakdown are well known, streamers are 
intrinsically the most difficult to investigate both experimentally and computationally, 
as they are very fast, very far removed from thermal equilibrium and sometimes even 
from equilibrium with the local electric field, and as they evolve on multiple scales in 
both space and time. The spatial scales range from the collisional mean free path of 
electrons, through the generic inner structure of a single streamer with its characteristic 
field enhancement, and on to the growth of interacting streamer trees consisting of 
tens-of-thousands of channels. The temporal scales range from electron impact collision 
times to the characteristic times for the formation of arcs, lightning and afterglow. An 
understanding of the full sequence of events in sparks or lightning must therefore include 
a detailed description of the streamer. Furthermore, the streamer by itself is active in 
corona reactors as well as in sprite discharges. Although various approximate schemes 
have been proposed for describing streamers, a detailed ab initio kinetic description at 
all spatial and temporal scales does not presently exist. 

1.2. Pertinent questions 

There arc a number of persistent theoretical questions concerning the nonlinear 
dynamical properties of streamers that have continued to elude fundamental description 
at the kinetic level, despite many decades of research. Among these are the necessary 
and sufficient conditions for their initiation, what determines their diameters and 
propagation speeds, when do streamers branch, how do they interact with one another 
or with dielectric walls, or with dust or water droplets or ice particles, when do they 
extinguish, how do streamers evolve across strong material or electrical gradients in the 
underlying medium, how do boundary conditions related to various kinds of electrodes 
affect the properties of streamers, whether these be solid laboratory electrodes or the 
ionosphere, how does the gas composition, its electronegativity and photoionizability 
and the level of background ionization, determine streamer structure and dynamics, and 
how are streamer properties affected by an anisotropic conductivity such as induced by 
a strong external magnetic field? 

All these questions have been grappled with in the past in one form or another by 
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the various research groups in their respective laboratory and geophysical investigations, 
but ad hoc approaches have often had to be used in the absence of a fundamental theory 
to provide key parameters. Examples of this approach are lumped circuit analogies to 
determine the current flow in the streamer body behind the head, or 1.5D computational 
models for streamers with pre-determined radius, or the dielectric breakdown model to 
describe the overall branched structure of many streamers. These have proven to be 
very useful concepts that embody considerable physical insight, but a complete streamer 
theory would be able to provide detailed descriptions of all these quantities, as well as all 
other quantities associated with a streamer, from first principles. It is the development 
of such a comprehensive theory that stands as one of the goals of streamer research. 

1.3. Methodological progress and the topics of the workshop 

In recent years, research methodology has progressed significantly as computational 
power and sophistication of algorithms has advanced, and as faster and more sensitive 
detection technology has become available: in plasma diagnostics and plasma modeling, 
in observations and modeling of atmospheric electricity, in electric pulsed power 
technology, and in multi-scale computing and nonlinear dynamics. The commonality 
and complexity of the breakdown problem is an incentive to cross disciplinary borders 
between plasma physics and chemistry, electrical engineering, industrial applications 
and geophysics, nonlinear dynamics and computational science. In fact, the frequency 
of cross-referencing among these disciplines in published research papers has increased 
considerably in recent years. 

The workshop at the Lorcntz Center in Leiden, The Netherlands, in October 2007 
created a forum for direct knowledge transfer among the various disciplines, and also 
for further development of methods and concepts within each discipline. It focused on 
the following topics relevant for laboratory, industrial and atmospheric discharges: 

(a) Microscopic mechanisms in different gases, their experimental identification and 
model implementation. This included photoionization, background ionization, electron 
impact excitation and ionization, thermal and nonthermal attachment and detachment 
processes, inelastic coUisions, etc. That the rates of the fast processes in the streamer 
head scale to good approximation with gas density, is a key to understanding why 
streamers at standard temperature and pressure and sprite discharges at high altitudes 
of the atmosphere are physically similar. 

(b) Macroscopic spatial structures. This includes observation and modeling of 
single streamers, their branching or extinction, and interactions among streamers, 
experimental methods for discharge observations in the laboratory and the atmosphere, 
and computational methods for studying their multi-scale nature. 

(c) The electron energy distribution at the streamer tip and inside the discharge and 
the subsequent plasma chemistry. Here, we include plasma chemistry in corona reactors 
and other industrial applications, the impacts of lightning and transient luminous events 
on atmospheric composition, spectroscopic signatures of ionization, and its potential for 
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X-ray emission in both the laboratory and in hghtning. 

(d) Electrical properties of the discharge and their electromagnetic emissions. 
This includes coupling to pulsed power circuits, measurements of inner electric fields, 
modeling of charge distribution and charge transport. 

(e) Processes on earlier or later time scales, on the one hand the initiation of 
streamers, and on the other hand later processes such as secondary streamers, transition 
to leaders, glows or arcs, and afterglow effects, heating and convection. 

2. The papers in the cluster issue and other recent literature 

The articles within the cluster issue are a representative sampling of recent research in 
various fields in which streamers play a role. We here briefly discuss how they relate 
to each other and indicate where they fit into the broader field of streamer research, 
without attempting a complete review. 

2.1. The structure of streamers in ambient air 

The first five articles in the cluster issue deal with the basic macroscopic characteristics 
of streamers in air at standard temperature and pressure; this includes their diameter, 
speed of propagation, and their branching into trees of many interacting streamers. 
The article [1] by Winands et al. presents measurements of many positive or negative 
streamers emerging together from a wire electrode, showing that streamers of different 
polarity largely resemble each other. The authors have reported similar observations 
earlier in p], [3]. In contrast, Eichwald et al. ^ and Nudnova et al. jS] present 
measurements and simulations of single positive streamers that emerge from a needle 
electrode; Eichwald for d.c. voltage and Nudnova for sharp voltage pulses similarly to 
the experiments of Winands. The focus on positive streamers is quite generic for the 
recent streamer literature. The measurements presented by Briels et al. [6] elucidate 
this apparent contradiction. They cover a voltage range of 5 to 95 kV with different 
power supplies and show that for low voltages only positive streamers emerge while for 
higher voltages, negative streamers are generated as well; with increasing voltage the 
negative streamers more and more resemble the positive ones; they both form multiply 
branched trees. 

The article of Nudnova et al. ^ builds on earlier work by Pancheshnyi and 
Starikoskii [71 El [9l [lO]; it focuses on simulation and experimental reconstruction of 
the head structure of positive streamers as revealed in optical emissions from the light- 
emitting shell around the streamer head. This inner head structure with its associated 
strong field enhancement at the streamer tip is the most important element of a streamer, 
and its properties play a pivotal role in the dynamics, as many simulations show. Its 
structure as determined in simulations is here confirmed in experimental observations. 

Within the voltage range of 5 to 96 kV, Briels et al. [6] find that the streamer 
diameter can vary by a factor of 15 and the propagation speed by a factor of 40; 
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they present detailed measurements of diameters and velocities as a function of the 
applied voltage, as well as an empirical fit of the velocity as a function of the diameter. 
Experiments are here characterized by voltages rather than by hypothetical average 
electric fields, as the fields in a wire-to-plane or needle-to-plane electrode geometry 
are quite inhomogeneous. The electric field near the needle or wire electrode where 
the streamers are launched is primarily determined by electrode geometry and applied 
voltage, and depends only weakly on the gap distance. This local field together with 
the availability of free electrons and with emission processes on the electrode surface 
determine the streamer inception; this inception process was recently imaged in 
After inception in the high field zone near the electrode, the streamers propagate into 
an undervolted region, i.e., into a region where the local electric field in the absence of 
streamers is below the breakdown value. 

The asymmetry between positive and negative streamers in air (that also plays 
in lightning physics [12]) is studied in the simulations by Luque et al. [13]. Here, 
streamers both in homogeneous and in inhomogeneous fields are studied. First double 
headed streamers in homogeneous overvoltage fields are investigated, similarly to earlier 
simulations by Liu et al. [TU |15]. Then positive and negative streamers are launched 
from the high field region near needle electrodes, and they propagate into a low field 
region as in the experiments [6]. Positive streamers propagate with a similar relation 
between velocity and radius as in the experiments [6], while negative streamers at low 
voltages broaden and extinguish, also in agreement with experiments. 

A characteristic of many streamers and sprites is that discharge channels appear in 
groups; they can branch or extinguish, and sometimes they reconnect. The full three- 
dimensional structure of streamer trees emerging from a needle electrode was recently 
analyzed through stereographic imaging, and the branching angles were determined by 
Nijdam et al. \T6\. Fully three-dimensional simulations of two interacting streamers 
were performed by Luque et al. [17] showing that equally charged streamer heads do 
not necessarily repel each other, but in air they also can merge due to the nonlocal 
photo-ionization reaction. 

2.2. The influence of magnetic fields and the role of medium and density 

Manders et al. [18j study short positive and negative streamers in ambient air in 
magnetic fields of 2.5 to 12.5 Tesla. They observe that the streamer path bends due to 
the Hall effect and they reconstruct the local electric fields. These measurements are 
also relevant for sprite discharges above 75-85 km where the thermal electron collision 
frequency drops below the cyclotron frequency in the geomagnetic field. 

Streamers are observed not just in air, but in almost any ionizable, nonconducting 
medium, be it gas, liquid or solid. We recall that phenomenological concepts of streamer 
propagation as presented in [19] were originally developed in the context of ionization 
fronts in semiconductors [20]. In the present cluster issue, the review by Kolb et 
al. [21] on streamers in water and other dielectric liquids is included. Without aiming 
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at any completeness of references, we also recall investigations in liquid nitrogen and 
helium [22| [23] and in gaseous nitrogen-oxygen mixtures [25| [26] . Discharges in 
other gases besides air are of interest; we refer to [271 EHl EH] for streamer processes in 
the ignition of high pressure discharge lamps, and to the collection of articles [30] and 
to [311 132] for a discussion on discharges in the atmospheres of other planets. Images of 
streamers in ambient air at 1 to 9 bar can be found in [33j . 

Similarity laws between phenomena in different gas densities such as streamers 
in ambient air and sprite discharges at high altitudes of the atmosphere are based 
on theoretical considerations [151 EH [35]. These similarity laws are reviewed and 
experimentally tested on positive streamers in [26] where the pressure was varied from 
0.013 to 1 bar and the gap distance from 10 to 160 mm; the extrapolation to sprites 
works very well. We remark that the relevant scales of electric fields, electron and ion 
densities, lengths and times depend on gas density, whereas the voltage does not [26l[35] . 
The paper by Stenbaek-Nielsen and McHarg [36] on high time resolution imaging of 
sprite discharges underscores the similarity of streamers and sprites from the geophysical 
side; the authors show the same light emission structure of the head, and analyze 
onset, velocities, brightness and size and branching, as well as later luminous and bead 
structures. 

The essence of the streamer mechanism is the self generated field enhancement at 
the streamer head. Adachi et al. [37] measure approximately the local field in sprite 
heads through a spectroscopic method. They estimate the charge moment transfer of 
the generating lightning stroke in the underlying thunderstorm, and deduce from it the 
evolution of the background electric field in the sprite region. This work further supports 
that indeed the streamer mechanism is active in sprite discharges. An important 
counterpart to the estimation of local electric fields from observations is the measurement 
of lightning currents and of the associated charge moment transfers between cloud and 
ground. Their global distributions and seasonal variations are discussed by Sato et 
al. [38]. 

The nonlocal photoionization reaction is an important ingredient for positive 
streamers in air specifically, but the actual reaction rates and photoabsorption lengths 
are not very reliably known, as discussed in [5l [T7]. (The different case of repetitive 
streamer discharges with memory effects is discussed by Pancheshnyi in [39]; the 
d.c. driven streamers studied by Eichwald et al. [1] are repetitive as well.) There 
are actually two macroscopic observations that can give an indirect handle on the 
role of photo-ionization or pre-ionization relative to impact ionization, namely the 
bending of the streamer in a strong magnetic field [H], and the merging or repulsion 
of two streamer heads [17]. The fact that photoionization is increasingly suppressed at 
pressures above approximately 80 mbar and for a decreasing oxygen-nitrogen ratio can 
be used in experiments like [6], [16], [251 [26] to study this interaction further on macroscopic 
phenomena. 
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2. 3. UV and X-ray emissions, chemical products and energy deposition 

A characteristic and important feature of streamer discharges and transient luminous 
events are UV and X-ray emissions, chemical products and energy deposition. 

X-ray emissions and gamma ray bursts from lightning are currently under very 
active investigation; the cluster issue contains two papers that investigate possible 
sources of intensive radiation. Nguyen et al. [IQ] measure multiple gamma ray 
bursts from a MV discharge in the laboratory. In contrast to earlier experimental 
investigations [HI they clearly attribute the bursts to the streamer- leader phase 
of the discharge. Milikh et al. propose that UV flashes observed from a satelhte 
are due to the streamer zone of upward propagating gigantic blue jets. Theoretical 
investigations [HI HSl HHl HTj HH] of this effect are reviewed in section 12. 4[ 

The total energy of sprites and other transient luminous events as inferred from 
three years of optical emission measurements made with the ISUAL instrument aboard 
the F0RM0SAT2 satellite are described Kuo et al. [49]. While the energy deposition 
per event is of the order of 50 MJ for sprites, halos and elves, most of the energy is 
deposited in the upper atmosphere by elves, as they are the most frequent. 

The streamer head can be considered as a self-organized, highly efficient plasma 
reactor that travels through the gas; it initiates various chemical reactions that are 
important both in technical applications and in atmospheric chemistry. Van Heesch 
et al. [50] report extremely efficient generation of O* radicals and ozone with rapidly 
pulsed streamer corona discharges when fast voltage pulses in the range of 50 to 90 kV are 
applied; they even reached a conversion of dissipated electrical energy into ozone above 
50%. O* radicals and ozone are widely used for disinfection purposes; its generation in 
other corona reactors is discussed in [2], HI [51] . 

Sprite chemistry is a new subject treated in the cluster issue by Gordillo- 
Vazquez [52lj following the seminal study of Sentman et al. [53]. These studies model the 
nonstationary, non-LTE chemical response of the upper atmosphere near 70 km to the 
passage of a streamer of typically observed magnitude. Using a chemical kinetics scheme 
involving several dozen neutral and ionic species and several hundred coupled chemical 
reactions, the models follow the production, loss and chemical transformation of species 
through a complex sequence of processes starting with breakdown in the streamer 
head and continuing until the system returns to predischarge equilibrium. The models 
include what are believed to be the principal processes involved in electrical discharges; 
ionization, excitation, dissociation, coUisional and radiative deactivation, attachment, 
detachment, recombination, charge exchange, and various chemical reactions. 

2.4. Simulations and theory 

Due to the many length and time scales involved and due to their nonlinear coupling 
in a dynamics far from equilibrium, streamer physics is a challenging multiscale 
problem [35] . The simulation of propagating single streamer in cylindrical symmetry in 
fluid approximation with its inner shell structure is after more than 20 years of study 
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now reaching maturity, a few recent papers are 01 El [71 [HI 13 [131 [IH [151 [311 [3S1 [Ml [S51 
[561 [571 [58] . The cluster issue only contains such studies, plus a study by Naidis [59] of 
repetitive discharges and the transition from streamer to the spark break-down; these 
long times can presently only be treated by analytical approximations. In the present 
section we also briefly review other recent developments in simulations and theory. 

An efficient computational treatment of the nonlocal photoionization reaction in 
oxygen-nitrogen mixtures like air for fluid models was recently proposed independently 
by Segur et al. [561 [SIl [60] and Luque et al. [58] ; the cluster issue contains a detailed 
follow-up study of different computational strategies for calculating the monochromatic 
radiative transfer in air in fluid models for streamers |61j . 

A vital computational issue is adaptive grid refinement that is capable of accurately 
resolving the inner spatial structure of the streamer head while avoiding the high 
computational costs of solving Laplace's equation in the large outer non-ionized region 
on the same fine grid as the in head region. Adaptive grid refinement in finite volume 
codes for single streamers was developed and thoroughly discussed by Montijn et 
al. [5l], it was extended by Luque et al. with photo- ionization [58] and to full three 
dimensions [17] . The cluster issue contains a study on adaptive grid refinement in finite 
element codes by Papadakis et al. [62j; another parallel adaptive mesh refinement code 
on Cartesian grids was recently presented by Pancheshnyi et al. [63]; Kolobov et al. also 
have recently published some streamer results based on dynamically adaptive Cartesian 
meshes [M]- Note that avalanches [55], propagating streamers and unstable branching 
streamers form a ladder of increasing computational complexity as density gradients 
become steeper and the physical branching instability should be clearly distinguished 
from a numerical instability. The branching time was found to saturate with decreasing 
size of the numerical grid in simulations by Montijn et al. [65]; these simulations were 
performed in cylindrical symmetry as many before and therefore in fact give an upper 
bound for the branching time in full three dimensions, but they do not follow the true 
evolution after branching. 

Full three-dimensional simulations of streamers in fluid approximation pose a major 
computational challenge. Interesting phenomena include the full branching structure 
after the streamer head has destabilized, the interaction of several streamers and almost 
all interactions of streamers with electrodes, walls, particles or local inhomogeneities. 
Kulikovsky [66] and Pancheshnyi [39] have presented images of 3D simulations of 
streamers at the moment of branching with low numerical accuracy on uniform grids. 
Luque et al. [T7] recently have developed a parallelizable 3D code with dynamically 
adapting grid refinement, and they have studied the interaction of two streamers with 
this method. 

The study of interfacial motion in nonlinear dynamics gives firm theoretical support 
for the statement that streamers can branch without any noise or disorder if the streamer 
head has developed into a dynamically unstable state. In Section 5 of ^35], this branching 
concept of a Laplacian instability is confronted with the old branching concept based on 
stochastic ionization avalanches that can be found in many textbooks and that can be 
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traced back to Raether in the 1930's. An extensive analysis of a moving boundary 
approximation for the deterministic fluid model in [671 EHl EHl [70] places streamer 
branching in the mathematical context of other branching problems in nature such 
as viscous fingering, dendritic solidification etc. 

Even though streamers can branch in fully deterministic fluid models, random 
motion of individual particles or a disordered background can trigger an inherent 
branching instability earlier than it would occur in a deterministic system. The influence 
of random initial conditions was studied in a ID model by Arrayas et al. [71]. The 
influence of "bubbles" of slightly different density on streamer branching in 2D was 
recently studied by Babaeva and Kushner [72], while they studied the influence of 
dust particles on streamer propagation in [7^. These studies are obviously relevant 
for discharges in liquids with bubbles [21] as well as in thunderclouds or in dust devils 
on Mars; these processes as well as sliding discharges along dielectric walls will not be 
further reviewed here. 

There are hundreds to tens-of-thousands of streamer heads in the laboratory [U |50] , 
or in the streamer zones of lightning leaders and of blue jets or in sprite discharges; to 
model them in fluid approximation is beyond the capabilities of the current generation 
of computers, therefore models must be reduced. Several simple electrodynamic 
approximations have been proposed, e.g., by Bazelyan and Raizer [19], and they are 
extended and used in [13]. In a search for such simple electrodynamic models, the heads 
of simulated single streamers are characterized by charge content, voltage, radius and 
velocity in [13]. On the other hand, the study [71] shows that a group of streamers 
can exhibit different electrodynamic behavior than a single streamer. In particular, the 
electric fleld in the streamer channels can be completely screened in a group of streamers, 
but not in a single streamer. 

On the other end of the hierarchy of length scales, the stochastics of single particles 
may play a role in the development and dynamics of streamers; ultimately, a fluid 
model is just an approximation of the true particle dynamics. In fact, a comparison 
of a Monte Carlo and a fluid realization of the same streamer ionization front shows 
that electron energies and the plasma density behind the front increasingly deviate 
with increasing electric fleld [45]. This is, of course, a pre-cursor of the electron run- 
away effect that could be a source of X-ray or gamma radiation of energetic streamers. 
This idea motivated Moss et al. [H], Li et al. [15] and Chanrion and Neubert [16] to 
investigate Monte Carlo models of streamers and sprites. Moss and Li did this in one- 
dimensional approximation while Chanrion studies three-dimensional streamers and uses 
so-called superparticles. However, Li et al. have demonstrated in another paper [47] that 
superparticle simulations of avalanches and streamers can suffer from strong numerical 
artifacts. They are developing a hybrid model [H] that treats the high fleld region with 
low electron numbers using a Monte Carlo model, and the ionized streamer interior with 
a fluid model. Monte Carlo models or hybrid Monte Carlo-fluid models can also track 
the primary excitations of molecules after the passage of the streamer front as an input 
for chemistry models, as well as trace fluctuations of particle numbers as a possible 
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additional trigger for the streamer branching instabihty. 
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